Boltzmann Test of Slonczewski's Theory of Spin Transfer Torque 
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We use a matrix Boltzmann equation formalism to test the accuracy of Slonczewski's theory of 
spin-transfer torque in thin-film heterostructures where a non-magnetic spacer layer separates two 
non-collinear ferromagnetic layers connected to non-magnetic leads. When applicable, the model 
predictions for the torque as a function of the angle between the two ferromagnets agree extremely 
well with the torques computed from a Boltzmann equation calculation. We focus on asymmetric 
structures (where the two ferromagnets and two leads are not identical) where the agreement pertains 
to a new analytic formula for the torque derived by us using Slonczewski's theory. In almost all 
cases, we can predict the correct value of the model parameters directly from the geometric and 
transport properties of the multilayer. For some asymmetric geometries, we predict a new mode of 
stable precession that does not occur for the symmetric case studied by Slonczewski. 



In 1996, Slonczewskii and Berger— predicted that an 
electric current flowing through a magnetic multilayer 
can exert a spin-transfer torque on the magnetic moments 
of the heterostructure. This torque can produce sta- 
ble magnetic precession and/or magnetic reversal, both 
of which have been widely studied experimentally and 
theoretically. 4 Figure ^ illustrates a common geometry 
known as a "spin-valve" where a non-magnetic spacer 
layer separates a thick "pinned" ferromagnetic layer from 
a thin "free" ferromagnetic layer. Non-magnetic leads 
connect the ferromagnets to electron reservoirs. 

Slonczewski^ developed a theory of spin-transfer 
torque that combines a density matrix description of the 
spacer layer with a circuit theory description of the re- 
mainder of the structure. He worked out the algebra 
for the case where the spacer layer is thin and the spin 
valve is symmetric (identical ferromagnets and leads) and 
found the torque L$ to be the same on the left and right 
spacer/ferromagnet interfaces. As a function of the angle 
between the two ferromagnets, 
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In this formula, / is the total current that flows through 
the structure, 
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i?l and Ri are effective resistances experienced by spin 
up and spin down electrons between the reservoir and the 
spacer layer. The conductance G = Ae 2 k 2 F /Air 2 Ti, where 
A is the cross-sectional area of the device. 

In this paper, we solve Slonczewski's equations for 
the general asymmetric case and derive formulae for the 
torques Ll(6) and Lr{9) on the left (x — xl) and right 
(x = xr) spacer/ferromagnet interfaces in Figure ^ We 
then compare these formulae with numerical results for 
the torque obtained from a matrix Boltzmann equation. 6 
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FIG. 1: Schematic of a five-layer spin-valve. A non-magnetic 
spacer layer separates two ferromagnetic layers whose mag- 
netizations are inclined from one another by an angle 9. A 
non-magnetic lead connects each ferromagnet to an electron 



reservoir. 



From this comparison, we are able to identify the physi- 
cal origin and systematic behavior of the effective resis- 
tances iZf and Ri in the interplay of inelastic scattering, 
spin-flip scattering, and interface scattering. For some 
asymmetric geometries, a previously unsuspected feature 
of the torque leads us to predict a new mode of stable 
precession that does not occur for the symmetric case. 

The electric current in a spin valve can be written as 
the sum of spatially varying currents carried by up and 
downi spin electrons: I = I-[(x) + Ii(x). The corre- 
sponding spin-current is Q(x) = I\{x) — The up 
and down spin voltage (chemical potential) drops along 
the spin valve are V\(x) and Vj,(x). Slonczewski^ con- 
tracted the function I^{x) to the values 1^ = I](xl) 
and 1^ — I-[(xr) and similarly for Ii(x). He also writes 
AVr = V-\(xr) — Vi(xr) for the up-down difference in 
the voltage drop from the right reservoir to a point in 
the spacer infinitesimally close to the interface between 
the spacer and the right ferromagnet. AVt, is defined 
similarly. These quantities are directly proportional to 
the spin accumulation at xr and Xl, respectively. 

With this model, Slonczewski wrote down (but did 
not completely solve) all the equations needed for the 
asymmetric geometry. In our notation, two of his linear 
equations relate the voltage drop differences to the spin 
currents: 







AV L (1 + cos 2 9) — G~ x Ql sin 2 - 2AV R cos 9 



= Ql (1 + cos 2 9) - GAVl sin 2 9 - 2Q R cos 9. (3) 

Two additional equations parameterize the voltage drop 
differences in terms of effective resistances Rl, Rr, Tl 
and tr: 

AV L = QlRl+Itl 

AVr = -QrRr-Itr. (4) 

Finally, Slonczewski derived expressions for the interfa- 
cial torques. At x = xr, the torque is 
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The torque Ll at the x — xl interface is (JHJ) with Qr 
and Ql exchanged. 

It is straightforward to solve the four equations © 
and I0J for the four unknowns, Ql, Qr, AVl, and AVr. 
From these, we find L R from JSJ to be 
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The parameters Pl, Pr, Xl and are defined in terms 
of R L , R R , r L and r R as P and A are defined in in 
terms of R and r. For the symmetric case, A^ = Xr, = A 
and Pl — Pr, = P. This makes q_ = and © reduces 
to Slonczewski's formula with Lr = Ll=Ls- We 
will see later that the term proportional to q_ in © can 
affect the magnetization dynamics of the spin valve in a 
qualitative way. 

To test we computed the torque on each 

spacer/ferromaget interface using a Boltzmann equation 
formalism^ The two ferromagnetic moments in Figure ^ 
are not collinear, so there is no natural spin quantiza- 
tion axis in the spacer layer. Therefore, we expand the 
deviation of the semi-classical electron occupation func- 
tion from its equilibrium value in a basis of Pauli spin 
matrices: 
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In the spacer layer, each of go, g x , g y , and g z satisfies 
a linear Boltzmann equation that takes account of the 
driving current and of resistive and spin-flip scattering 
in each material layer. In each ferromagnet and in the 
adjacent lead, it is sufficient to use go and g z referenced to 
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FIG. 2: Spin-transfer torque at x = xr for a spin-valve with 
layer thicknesses 5 nm/40 nm/1 nm/1 nm/180 nm. Solid cir- 
cles are Boltzmann equation results. Solid curve is Eq. HJ de- 
rived from Slonczewski's theory. Dashed curve is Eq. HJwith 
all bulk scattering removed. 

the fixed direction of magnetization in each. We assume 
a spherical Fermi surface for each material. We adopt 
a one-dimensional approximation g(k, r) = g(k, x) and 
stitch together the solutions in each layer with match- 
ing conditions that reflect the complex, spin- and wave 
vector-dependent reflection and transmission amplitudes 
of each interface. The latter are determined from a pre- 
viously published parameterization, 9 

Electrons that enter a non-magnetic lead from the ad- 
jacent reservoir are assumed to have an equilibrium occu- 
pation function. For our numerical work, we have chosen 
material and geometry parameters typical of experiments 
performed on Cu/Co/Cu/Co/Cu spin valves. The spe- 
cific numerical values used can be found in Ref. 6. With 
the final occupation function in hand, it is straightfor- 
ward to compute the spatially varying voltage, spin accu- 
mulation, spin current, and spin-transfer torque. Details 
will be published elsewhere^ 

The filled circles in Figure |2] show a typical Boltzmann 
result for the spin-transfer torque at x = xr as a function 
of the angle 9 between the two ferromagnets for an asym- 
metric geometry. The solid curve in Figure [21 is the same 
quantity, Lr(9), computed from JHJ. The agreement is 
excellent, as it is for essentially all other geometries we 
have studied with thin spacer layers. For comparison, 
the dashed curve is the Slonczewski torque JBJl with the 
bulk scattering removed. This "interface-only" situation 
is manifestly symmetric (g_ = 0). 

For the solid curve plotted in Figure^ the relative im- 
portance of the two terms in JHJ is q-/q+ — 0.27. We 
find that this ratio does not exceed about 0.5 for physi- 
cally sensible geometries. We will address the qualitative 
consequences of g_ ^ at the end. First, we describe 
our method to determine the torque parameters A^, A^, 
Pr, and Pl that produced the solid curve in Figure |21 
We begin by writing an exact expression for the voltage 
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I a (x) is I a (x) minus the corresponding current in an in- 
finite bulk sample of the material present at point x. 
The average resistivity in |(?JJ| is p — (pj + p\)/2. We 
will also need the resistivity difference Ap = (pi — p{)/2. 
Both p and Ap contain delta functions at the four non- 
magnet/ferromagnet interfaces to take account of spin- 
dependent interface scattering. 

To make progress, we use a drift-diffusion resultii to 
the effect that /-[ (x) and jj_ (x) decay exponentially close 
to the x = xq interface (Figure ^) in both directions. 
The decay length is the spin-flip length in each mate- 
rial, whether ferromagnet (F) or non-magnet (N). In that 
case, an approximate expression for JHJ is 

AV L = QoPn 4 + Q Q p F 4 + /a Pf dp 

+ AVi + AVc. (10) 

In ltTU|) . Q = Q(xo), AVi and AVc are voltage drops at 
the internal interfaces and at the reservoir contact, and 

4 = lf { [1 - cxp(-iM)] 

4 = ® [l-exp(-tM)]. (11) 

The effective lengths (|ll|l appear because, due to spin- 
flip scattering, only electrons within dp or d^ of the 
ferromagnetic interfaces can accommodate the dissimi- 
lar spin-currents characteristic of the ferromagnets and 
the non-magnets in equilibrium. 

The relationship between Q and Ql is non-trivialii 
except when the ferromagnet is very thin (£p <C iff). In 
that case, Qo — Ql, and we can connect (|1U|) to an d 
© to get 

A| = G(p N d^ + p F 4 + Ri + R c ) 

P L = Al 2 (Ap F 4 + AR!). (12) 

These two formulae (and similar ones for A# and Pr), 
together with © and Q are the principal results of this 
paper. We used to compute the solid curve in Fig- 
ure El The interface resistance GRi ks 0.97 and contact 
resistance GRc ~ 1.1 were extracted from the Boltz- 
mann solution. GRi differs from the experimental value 
by about 15%i£ 

The parametrization 10 is well suited to study the 
behavior of L{9) when we vary the geometry of the mag- 
netic heterostructure and the material parameters in our 
Boltzmann calculations. For convenience, we did this for 
symmetric geometries. Figure |21 confirms that A 2 is a 
linear function of £n when ifl 3> £n but saturates when 
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FIG. 3: Torque parameter A 2 as a function of £n for large and 
small values of lf f . 
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FIG. 4: Torque parameter A 2 as a function of tn for different 
values of lf t . In all cases, 1^ = oo. The curve for lf f S> tp 
corresponds to A 2 oc t$t. 

% ~ Igf- Interestingly, the saturated value of A 2 varies 
linearly with lf f — l N (Z N is the inelastic scattering length) 
rather than with lf { as predicted by ifTTJl . For long leads, 
this can be understood from the fact that conventional 
resistive scattering is needed to build up non-equilibrium 
spin accumulation in the non-magnet while spin-flip scat- 
tering works to return the non-magnet to equilibrium. 

Figure0]shows the variation of A 2 with lead length for 
different values of the spin-flip length in the ferromagnet. 
This calculation puts lf { — * oo, so we expect from Ijllf) 
that A 2 oc £n- This is indeed the case when lf f 3> t F . 
However, when the iff is comparable (or less than) the 
ferromagnetic layer thickness, the torque parameter sat- 
urates. This is a signal that our approximation Q ~ Ql 
has broken down. In this limit, fast spin- flipping in the 
ferromagnet reduces Qq to a value much less than Ql- 
When that is the case, relatively little spin-flip volume in 
the non-magnet is needed to reduce the spin current in 
the non-magnet to zero (its equilibrium value). 

A characteristic difference between Ls(0) in Q for a 
symmetric geometry and Lr(9) in © for an asymmet- 
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FIG. 5: Left: Spin transfer torque and damping torque for a 1 nm/40 nm/1 nm/1 nm/1000 nm spin valve (q-/q+ ~ 0.36); 
Middle: the total torque as the current increases; Right: the angle between two ferromagnetic moments as a function of current 
I in arbitrary units. 



ric geometry can be seen from the difference between the 
dashed curve and the solid curve in Figure [21 The for- 
mer has a bump (maximum) in the interval it/2 < 8 < tt 
only. The latter has a small additional bump in the in- 
terval < 8 < ir/2 that comes from the <?_ termi^ This 
small change is enough to produce stable magnetization 
precession for some asymmetric geometries. 

Consider a spin valve in the presence of an external 
magnetic field aligned with the magnetization of the thick 
ferromagnet. If we ignore shape anisotropy and lattice 
anisotropy, the total torque acting on the thin ferro- 
magnetic film when electrons flow from right to left in 
Figure ^ is the sum of the spin-transfer torque Lr(8) 
and a Gilbert damping torque "fH sin 8 (left panel of Fig- 
urc[5J). As the current increases, the spin-transfer torque 
increases and eventually destablizes an initial state with 



parallel moments. Stable precession occurs at angles 
where the total torque changes from positive to negative 
(middle panel of Figure EJ). When the total torque be- 
comes everywhere positive, the system abruptly switches 
to the anti-parallel configuration (right panel of Fig- 
ure ISJ) • There is no regime of stable precession if the 
zero-current state is anti-parallel. 

In summary, we have shown that Slonczewski's theory 
of spin-transfer torque in spin-valves can reproduce 
the results of Boltzmann equation calculations when 
the non-magnetic spacer layer is thin. When the 
ferromagnetic layers are also thin, the parameters 
of the theory can be calculated from first princi- 
ples. One of us (J.X.) acknowledges support from the 
National Science Foundation under grant DMR-9820230. 
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